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FRACTIONAL TERM STRUCTURE MODELS: 
NO-ARBITRAGE AND CONSISTENCY 

By Alberto Ohashi 1 

Ibmec-Sdo Paulo and Universidade Estadual de Campinas 

In this work we introduce Heath-Jarrow-Morton (HJM) interest 
rate models driven by fractional Brownian motions. By using support 
arguments we prove that the resulting model is arbitrage free under 
proportional transaction costs in the same spirit of Guasoni [Math. 
Finance 16 (2006) 569-582]. In particular, we obtain a drift condition 
which is similar in nature to the classical HJM no-arbitrage drift 
restriction. 

The second part of this paper deals with consistency problems 
related to the fractional HJM dynamics. We give a fairly complete 
characterization of finite-dimensional invariant manifolds for HJM 
models with fractional Brownian motion by means of Nagumo-type 
conditions. As an application, we investigate consistency of Nelson- 
Siegel family with respect to Ho-Lee and Hull- White models. It turns 
out that similar to the Brownian case such a family does not go well 
with the fractional HJM dynamics with deterministic volatility. In 
fact, there is no nontrivial fractional interest rate model consistent 
with the Nelson-Siegel family. 

1. Introduction. Financial models driven by semimartingales and Markov 
noises have been intensively studied over the last years by many authors. In 
general, absence of arbitrage is the basic equilibrium condition which fulfills 
the minimum requirement for any sensible pricing model. On the other hand, 
empirical studies propose models which are not consistent with this basic 
assumption. In particular, some evidence of nontrivial long-memory behav- 
ior in bond markets has been recently suggested by many authors (see, e.g., 
[21] and other references therein). In most cases, the presence of long-range 
dependence in short-rate interest rates seems to be common and it is origi- 
nated by the fundamentals of the economy. In this regard, it is important to 
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study bond markets with extrinsic memory driven by non-Markovian noises 
which allow nontrivial long-range dependence over time. 

Recall that in the classical Musiela parametrization the forward rate rt 
satisfies a stochastic partial differential equation (henceforth abbreviated by 
SPDE) of the following form 

(1.1) dr t (x)=^r t (x)+a H m(t,rt(x))^dt + Y^v j (t,r t ( y x))dB{, 

where anjM is the so-called Heath-Jarrow-Morton (henceforth abbreviated 
by HJM) drift condition which is completely determined by the volatilities 
(o~ J )j>i under a risk-neutral measure, (B J )j>i is a sequence of stochastic 
noises and x is the time to maturity. The forward rate rt is considered as 
a Hilbert space-valued stochastic process. Due to this infinite-dimensional 
intrinsic nature, it is important to understand the relation between forward 
curves x i— > rt (x) at time t > and finite-dimensional parameterized families 
of smooth forward curves, frequently used in estimating the term structure 
of interest rates (e.g., Nelson-Siegel and Svensson families). 

Originally proposed by Bjork and Christensen [3] and recently studied by 
Filipovic and Teichmann [10, 11], the so-called consistency problems refers to 
the characterization and existence of finite-dimensional invariant manifolds 
with respect to t \— > rt- In fact, the stochastic invariance is essentially equiv- 
alent to a deterministic tangency condition on the coefficients J| , ohjm j o~ 
in (1.1). In particular, if the short-rate of interest rates exhibit long-range 
dependence then standard statistical procedures may be misspecified, since 
in this case the classical HJM no-arbitrage drift restriction may not be the 
correct one. Moreover, by fixing (nonsemimartingale) long-memory stochas- 
tic noises (Bj)j>i in (1.1), one has to obtain new tangency conditions on 
the coefficients of (1.1) to get appropriate arbitrage- free invariant parame- 
terized families of smooth forward curves. This is the program that we start 
to carry out in this work. 

We have chosen the driving noise in (1.1) given by the fractional Brownian 
motion (henceforth abbreviated by fBm) with Hurst parameter H € (1/2, 1). 
For many reasons (see, e.g., [26]), the fBm appears naturally as the canonical 
process with nontrivial time correlations inserting memory into system under 
consideration. Indeed, the main difficulty in dealing with fBm is the fact 
that such process is a semimartingale if and only if it is a standard Brownian 
motion (H = 1/2). This lack of semimartingale property immediately implies 
that fBm allows arbitrage opportunities (for any H 7^ 1/2) in the absence 
of transaction costs. In particular, Gapeev [12] obtains an explicit arbitrage 
strategy via pathwise integrals in a frictionless bond market driven by a fBm 
with H > 1/2. 

The main goal of this work is to introduce arbitrage-free HJM interest 
rate models driven by a fBm under arbitrary small proportional transaction 
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costs in the bond market. In this paper, the forward rate is considered as 
the solution of a SPDE (in Skorohod sense) of type (1.1) under the Musiela 
parametrization. In this work, we only treat the case of deterministic volatil- 
ities, leaving open the general stochastic volatility case for future research. 

Under deterministic volatility assumption, we obtain a drift condition 
which is similar in nature to the classical HJM no-arbitrage drift restric- 
tion. Although such condition is not sufficient to ensure no-arbitrage in the 
market, when combined with an additional mild condition on the volatilities 
it results in absence of arbitrage in the same spirit of the works [14, 15], 
where the support of the driving noise plays a key rule in the no-arbitrage 
characterization for markets with transactions costs. 

In the second part of this paper, we characterize finite-dimensional invari- 
ant submanifolds for HJM models driven by fBm by means of Nagumo-type 
conditions. Such characterization is the key ingredient to tackle the consis- 
tency problems related to the model. As an application of these abstract 
results, we investigate consistency of the Nelson-Siegel family with respect 
to Ho-Lee and Hull- White models driven by fBm. Similar to the Brownian 
case, such a family is not consistent with respect to these models. In general, 
we arrive at the same classical result of the Brownian case: no nontrivial in- 
terest rate model with deterministic volatility structure is consistent with 
Nelson-Siegel family. 

This work is organized as follows. In Section 2, we define the basic SPDE 
which describes the term structure of the interest rates. In Section 3, we 
give some general results regarding portfolios and absence of arbitrage in 
fractional bond markets. In Section 4, we characterize finite-dimensional in- 
variant forward manifolds with respect to HJM models driven by fBm. In 
Section 5, we examine consistency of the Nelson-Siegel family with respect 
to concrete interest rate models. The Appendix contains some technical re- 
sults with respect to the integrals which appear in the self-financing trading 
strategies introduced in Section 3. 

2. Fractional term structure equations. Throughout this paper we are 
given a d-dimensional fBm (3 = (/3- J )^ =1 with parameter 1/2 < H < 1 on 
a stochastic basis (J7, (J r j)t>o,P) satisfying the usual conditions. In other 
words, it is a centered <i-dimensional Gaussian process with continuous sam- 
ple paths, [3o = and covariance 

E(Pi-[3lM-%) = 5 l3 \t- s \ 2H ~ 2 

for t, s 6 M+ and i, j = 1, . . . ,d. For a detailed discussion on the stochastic 
analysis of the fBm, the reader may refer to [18, 26]. 

We denote by Ct the space of real- valued continuous functions on a metric 
space T. In what follows, we consider the following subset of M 2 

A 2 := {(t,T) G M 2 |0 < t < T < oo}. 
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Let us consider a term structure of bond prices {P(t,T);(t,T) G A 2 } 
where P(t, T) is the price of a zero coupon bond at time t maturing at 
time T. We assume the usual normalization condition 

P(t,t) = l Vi>0, 

and P(t,T) is a.s. continuously differentiable in the variable T. In this way, 
we introduce the term structure of interest rates {f(t,T); (t,T) G A 2 } given 
by the following relation: 

(2.1) P(t,T)=exp(- jT f(t,u)du\, (t,T) G A 2 , 

for some a.s locally integrable function / : [t, +00) — > M, the time t forward 
curve. In this paper, we adopt the Heath-Jarrow-Morton framework [17] 
in the fBm setting in terms of the Musiela parametrization where rt(x) := 
f(t,t + x) for (t, x) G and x = T — t is the time to maturity. 

In particular, we seek the forward curve x *— > rt(x) as a Hilbert space- 
valued stochastic process described by a linear SPDE 

d 

(2.2) dr t = (Ar t + a t )dt + J2<dPi r (-)=^E 

i=l 

in a separable Hilbert space E to be defined. The first-order derivative opera- 
tor A := ■4- is the infinitesimal generator of the right-shift family of operators 
{S(t);t > 0} acting on E. Of course, the SPDE (2.2) must be interpreted in 
the integral form. Moreover, in this work the stochastic integral is considered 
in the Skorohod sense [9] as a Paley-Wiener integral. 



2.1. The specification of the model. In this section, we specify a term 
structure model driven by a d-dimensional fBm. Let us assume for the mo- 
ment that the forward rate is given by the following system of stochastic 
differential equations 



( 



2.3) f(t,T) = f(0,T)+ fais^ds + Y, f <? l (s,T)df? s . 

Jo f^Jo 



From now on the coefficients and a are deterministic func- 

tions. Equation (2.3) is well defined if for each i = 1, . . . , d 

r T i-T t-T 



f \a(s,T)\ds+ [ I \o- i {s,T)\\<j i (t,T)\(l> H (t-s)dsdt< 00 
Jo Jo Jo 

for all < T < 00, where <j> H {u) := H(2H - l)\u\ 2H ~ 2 , ueR. 
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Let {S(t);t > 0} be the semigroup of right shifts defined by S(t)g(x) := 
g(t + x) for any function g:R + -> R. Fix (t,x) € R+. Then (2.3) can be 
written as 

f(t,t + x)=S(t)f(0,x) + [ S(t- s)a(s,s + x)ds 

(2.4) 

d ft 

+ Y. / S{t-s)a l {s,s + x)dPl. 
i=i Jo 

In (2.4) we deal with the parametrization T = t + x. The operator S(t) 
acts on f(0,x), a(s,x + s) and a 3 (s,x + s) as functions of x. By setting 

r t{x) :=f(t,t + x), 

it follows that 

P(t,T)=exp|-^ T *r t (x)da;|; (t,T)eA 2 . 

We can work out in an axiomatic way the minimal requirements on a 
Hilbert space E such that (2.4) can be given a meaning when 

(2.5) r t (-) = /(*,*+■); i£R+ 

is considered as the mild solution of (2.2). The strategy in finding a suitable 
state space E follows very similar to Filipovic [10]; see Section 4.2 and 
Theorem 5.1.1. So we omit the details and the reader may refer to this 
work. 



Remark 2.1. Recall that the minimal requirements on the state space 
E are the following ones: the right-shift semigroup is a Co-semigroup on E 
and the evaluation mapping (5 x ,h) = h(x) is a bounded linear functional on 
E. We then choose the state space E as defined in [10], Section 5. 

In the remainder of this paper, U is a (/-dimensional vector space with 
an orthonormal basis (ej)f =1 and C(U, E) is the space of bounded linear 
operators from U into E with the usual norm || • ||. We make use of the 
following notation: we set at(-) := a(t,t + •) and a = (<r J )y =1 , where 

4(-) :=o-tej(.) :=<t j (M + -); j = l,...,d. 

In this paper, we are interested in Gaussian interest rate models where 
we assume that the coefficients a : R+ — > E and a : R+ — > C(U, E) satisfy the 
following set of assumptions: 

(2.6) / ||a s ||£ ds + / 1 1 <r s 1 1 2 < oo for every < T < oo. 
Jo Jo 
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To ensure existence of a continuous version for the mild solution of (2.2) 
we assume there exists 7 G (0, 1/2) such that 

rT r T 

v~"' \\S (u)a u \\ \\S(v)(t v \\4>h(u — v) dudv < 00 

(2.7) 

for every < T < 00. 



In order to get a well-defined expression for the bond prices {P(t, T); (t, T) € 
A 2 } we also assume the following growth conditions: 



(21 



(2.9) 



/ ll cr tt( s )llR<i||cr t ,(r)|| Rt i^_H-(M — v)dudvdsdr < 00 
J[0,T]* 

for every < T < 00; 
Wu(t)\\Rd\\a v (t)\\ R d(j)H(u ~ v) dv dudt < 00 

for every < T < 00. 



[0,T] 



One should note that (2.6) yields J T ||5(t)<rt || 2 dt < 00 for every < T < 
00, and therefore the stochastic convolution is a well-defined E- valued Gaus- 
sian process given by 

£ f S(t-s)aidf3i, t>0. 

We assume the existence of a traded asset that pays interest. In other 
words, the unit of money invested at time zero in this asset gives at time t 
the amount 



S (t) : =exp|^ r s (0)ete}, 



where rj(0) = f(t,t) for t > 0. By considering Sq as a numeraire, the dis- 
counted prices are then expressed by 

(2.10) Zt (T):=^^-, (t,T) £ A 2 . 



In the sequel, we make use of the following notation: if v : [0, T] x 



is locally integrable in M + then Z v (s,T) := Jq s v s (x)dx. Under the above 
assumptions one can easily show the following lemma. 

Lemma 2.1. Assume that the coefficients a and a satisfy the assump- 
tions (2.6), (2.7), (2.8) and (2.9). Then the forward rate rt is the continuous 
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mild solution of (2.2). Moreover, the term structure of bond prices is given 
by the continuous process 



Proof. Equation (2.11) is a straightforward application of stochastic 
Fubini theorem (see [20]) in the fBm setting by using conditions (2.8) and 
(2.9). Conditions (2.6) and (2.7) allow the existence of a continuous mild 
solution of (2.2) as in [9]. □ 

3. The bond market: Portfolios and no-arbitrage. In this section, we 
discuss the basic setting for a fractional term structure of bond prices with 
transaction costs. It is well known that not every choice of stochastic integral 
makes sense in Finance (see, e.g., [5]). Therefore, at this point it is necessary 
to discuss if our choice of stochastic integral presents conceptual problems. 
We recall that integrals with respect to Cr + -valued processes play a key rule 
in the bond market theory. The integration theory for Banach space-valued 
semimartingales can be developed in many ways (see, e.g., [4, 6, 8]). In 
general, the fundamental idea is to make use of duality arguments between 
a "state variable" St (which represents a discounted price curve at time t in 
Ct) and a dual C^-valued process fi which represents a portfolio strategy. 
We stress here that this procedure is fully related to the intrinsic mechanism 
of a bond market and not to a given stochastic dynamics on S. In this way, 
we may adopt the same procedure in the nonsemimartingale case. 

Since we are dealing with a continuous Gaussian HJM model we shall 
make use of a mild integrability assumption [see (HI)] to define the cor- 
respondent wealth processes without any use of the fine structure of the 
Skorohod integrals which drive Z. In particular, by a standard limit proce- 
dure we avoid the Wick products in the definition of self-financing portfolios 
(see, e.g., [5] for a discussion on stock markets). For us, the most important 
property is an integration by parts formula (see Proposition A.l) in the same 
spirit of the works [14, 15]. 

We now introduce the notions of admissible self-financing portfolios in 
our context. Let us denote by Mt* the space of (finite) signed measures on 
[0, T*] endowed with the total variation norm || • ||tv- Let \i be a measure- 
valued elementary process of the form 



P(t,T) = P(0,T) 



(2.11) 




for (t,T) <E A 2 . 



N-l 



(3.1) 



Mt(w,-) := XF l x(t i ,t i+1 ](w,t)m i (-) 



i=0 
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where m ; <G Mr* , = t < ■ ■ ■ < T N < T* and F t G T u . 

We denote by Sb the set of elementary processes of the form (3.1) endowed 
with the following norm: 

(3.2) ||/z||y:=E sup Hxv- 

0<t<T* 

From now on all economic activity will be assumed to take place on the 
bounded set [0,T*] 2 . So we assume that Z t (T) = if (t,T) £ [0,T*] 2 . Under 
the hypotheses (2.6), (2.7), (2.8) and (2.9), the discounted price process 
Zt(T) satisfies the following condition: 

(HI) {Zt(T)\ (t,T) £ [0,T*] 2 } is a jointly continuous real-valued stochastic 
process such that 

E sup |^(T)| 2 <oo. 
(*,T)e[0,T*] 2 

If \i € Sb is given by (3.1) then we define 

rt N-l 

(3.3) / fj, s dZ s := V XFi(Z ti+lA t ~ z tiAt)mi, 

where Zt l mt i is the usual dual action. By Holder's inequality it follows that 



(3.4) E sup 

0<t<T* 



rt 

/ fi s dZ s < \\n\\ v E 1/2 sup \\Z S - ZtW 2 ^ < oo, 

JO 0<s,t<T* 



where || • ||oo denotes the usual (uniform topology) norm on the space of 
real-valued bounded functions defined on [0,T*]. Let V be the completion 
of Sb with respect to (3.2). By the estimate (3.4) and the definition of V we 
may easily define J fj, s dZ s for every [lEV. 

In the sequel, we denote Vt* the set of all partitions of [0, T*]. We also 
need the following assumption: 

(H2) U T *(fi) :=sup ne7 , T „ Euen \\Vk +1 ~^||tv is square integrable. 

By taking into account proportional transaction costs in the bond market, 
the liquidation value of a portfolio with zero initial capital is 

Vtiv) ■= zZ XF t { z U +1 ht ~ Z UM )mi 
U<t 

~ k Y z uW +1 -kZ t \tH\, 

where k is an arbitrary positive number and | • | denotes the total variation 
measure. The first term accounts for the capital gain of holding an elemen- 
tary strategy /i of the form (3.1) (without transaction costs) during the 
interval [0,t\. The second and third term account for the transaction costs 
incurred in various transactions and the eventual liquidation value of the 
portfolio, respectively. 
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By passing from a finite number of transactions to continuous trading 
one can easily show that if \x E V satisfies assumption (H2) then the above 
quantities converge to the following: 

(3.5) V t k (fi):= [* fi s dZ s -k f Z a d\n 8 \-kZt\iH\- 

Jo Jo 

In other words, under assumptions (HI) and (H2) it follows that any wealth 
process is the limit of elementary strategies in the sense that 

lim E sup \V t k {fi n ) -V t k (p)\ =0; k > 0, 

n ^°° 0<t<T* 

where /i n is a sequence of elementary strategies such that fi n — > [i in V as 
n — > oo. See the Appendix for more details, including the definition of the 
second integral in (3.5). 

Remark 3.1. One may also introduce another approach to define the 
integral for measure-valued integrands; the integration theory is reduced 
via Fubini theorems where the stochastic differential is interpreted in the 
Stratonovich sense [24]. In this case, J fi s dX s is replaced by 




where dX t (T) = &(T) dt + £i=i p\{T) is a SDE satisfying (HI). One 
can easily show that under additional standard Holder assumptions on the 
paths of n € V and a (see, e.g., [2, 24]), we do have Jq \i s dZ s = (/io Z) t for 
< t < T. 

Now we are able to introduce the following notions: 

Definition 3.1. We say that /x G V is an admissible trading strategy 
if it satisfies (H2), it is (weakly) JT t -adapted and there exists a constant 
M > such that V t k (n) > —M a.s. for every t <T*. An admissible trading 
strategy is an arbitrage opportunity with transaction costs k > on [0,T*], 
if V^*(/i) > a.s. and P{Vj?*(^) > 0} > 0. Therefore, the bond market is 
k-arbitrage free on [0,T*] with transaction costs k if for every admissible 
strategy \x, Vj?* (fi) > a.s. only if Vj?* (fi) = a.s. 

Remark 3.2. Since the main dynamics takes place on A 2 we do as- 
sume that the support of fi s is concentrated on [s, +co) for every admissible 
strategy fi. 

It is straightforward to prove the following result in the same spirit of 
Guasoni ([14], Proposition 2.1) by using the integration by parts formula 
(A.l). So we omit the details. 
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Proposition 3.1. Let us fix k > 0. If for every (Ft) t>o~ stopping time r 
such that P{t < T*} > we have 

Z T (r) 



(3.6) P<{ sup 

.r<t<T<T* 



1 



<fc,r<r >o, 



Zt(T) 

then the bond market is arbitrage free on [0, T*] with transaction costs k. 

One can also show by using similar arguments from Guasoni, Rasonyi 
and Schachermayer [15] that the fc-arbitrage-free property in Definition 3.1 
is essentially equivalent to the existence of a /c-consistent price system. Thus 

(3.6) is also a sufficient condition for it. In fact, a sufficient condition for no- 
arbitrage is the conditional full support property for Z which is equivalent 
to full support only if Z is Markovian. Hence in the non-Markovian setting 
it is more natural to find mild conditions on the volatility in such way that 
logZ has only full support. See Lemmas 3.1 and 3.2 for details. 

3.1. Absence of arbitrage. We prove that under suitable conditions on 
the volatility a = (a^)j =1 , the bond market model is /c-arbitrage free for 
every k > 0. The main ingredient in the no-arbitrage argument consists in 
the full support property on C A 2 , where A^» := {(t,T);0 <t<T< T*}. 
This property together with a suitable choice on the drift will result in 
fc-no-arbitrage for every k > 0. Recall that if X is a Polish space then a 
random element £ : $7 — > X has P-full support when P^ := P o (U) > for 
every nonempty open set U in X . 

Lemma 3.1. Let Y: Vl — > C A 2 be a measurable map such that X:=logY 

T* 

has F-full support. Then Y satisfies assumption in Proposition 3.1. 

Proof. Given e > and r a ^-stopping time such that P{r < T*} > 0, 
it is sufficient to check that 

r< sup \X(t,T)-X(T,r)\ <e,r<T* >0. 

l r<t<T<T* ' 

If p G then triangle inequality yields 

{ sup \X(t,T)-p(t,T)\<e/2,T<T*\ 

C sup |X(i,T)-X(r,r)|< 

l r<t<T<T* ' 

Let us consider V the set of polynomials p on A^* with rational coefficients 
such that p(0,0) = 0. We claim that there exists p^V such that 

(3.7) P{ sup \X(t,T)-p(t,T)\<e/2,T<T*}>0. 
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Suppose that (3.7) is violated for every p € V. Then we obtain 

{ sup \X(t,T)-p(t,T)\<e/2,T<T*\ 
(t,T)eAj„ 

C{r>T*} P-a.s.Vpe'P. 

Therefore 

(3.8) |J{ SU P I X (*» T )-P(*> T )I < e / 2 } C{r>T*} P-a.s. 

By the density of V in C A 2 ^ and the full support of X it follows that 

p(|J{ sup |X(t,T)-p(i,r)|<e/2}} = l 
l peP L (i,r) G A2, J J 

and therefore P{r < T*} = which is a contradiction. □ 

Remark 3.3. Recall that the fBm has 7-Holder continuous paths a.s. 
for any 7 < H. Moreover, one can prove the existence of the fBm Wiener 
measure on a separable Banach space W continuously imbedded on the space 
Ck + such that the elements of W are 7-Holder continuous functions on any 
compact interval. See [16] for the proof of this fact. 

The following remark turns out to be very useful for the approach taken 
in this work. 

Lemma 3.2. Assume that I a j(t,T) is X-Holder continuous on for 
every j > 1 where 1/2 < A < 1. Then the process Y^ d j=i Jo -^o-j ( s ) T) has 
F-full support on C^2 . 

PROOF. Fix (£ J )j =1 a sequence of 7-Holder continuous functions on 
[0,T*] where 1/2 < 7 < H. We recall that if l aj (t,T) is A-H61der continuous 
on Ay* then the pathwise Young integral J I a j(s,T) <i£| is well defined and 
there exists a constant C > which depends only on T* , 7 and A such that 



(3.9) 



Jo 



where || • \\ v denotes the usual 77-Holder norm. 

Moreover, the pathwise Young integral coincides with the symmetric inte- 
gral in Russo and Vallois [24] . Recall that we are assuming that the volatil- 
ities are deterministic functions and therefore the Gross-Sobolev derivative 



12 



A. OHASHI 



of X a j(t,T) vanishes for each j > 1 and (i,T) £ Ay*. Since the fBm has 7- 
Holder continuous paths a.s., Proposition 3 in [2] implies that the Skorohod 
integral 

fleiis^dPi; j>l, 
Jo 

can be interpreted as a pathwise Young integral. By the estimate (3.9) and 
Remark 3.3 it follows that each J l a j(s,-)dPi has F-full support on C A 2^. 
Moreover, since (P J )j>i is a sequence of real- valued independent fBm we 
then conclude that 

(*,T)-X; flAs,T)dP s 

j=l J ° 

has P-full support as well. □ 



By Lemma 3.1 and Proposition 3.1 we know that if logZ has P-full sup- 
port then the bond market is fc-arbitrage free for every k > 0. Moreover, if 
the volatility a = (a^)j =1 satisfies the assumptions in Lemma 3.2, there are 
infinitely many choices of a which give the full support property for logZ 
and therefore absence of arbitrage in the fractional bond market. But there 
is a canonical choice for the drift which gives the desirable property (see 
Definition 3.2, Remark 3.5 and Proposition 3.2): for every < T < 00 

(3.10) EZ t (T) =P(0,T), 0<t<T. 

As a direct consequence of Lemma 2.1 we have the following basic result. 



Corollary 3.1. Condition (3.10) holds if and only if the drift a sat- 
isfies the following equality 

i=i 1 h 

(3.11) 

+ ^ ai(y)dyJ o ^(- + t - 0)<f> H (t - 9) del. 
Proof. By Lemma 2.1 it follows that 

Zt (T) = P(0, T) exp j - j\ a ( 8 , T) ds - £ jT* l ai (s, T) d(i{ | ; 

(t, T) G A 2 . 
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If < T < oo, then the Ito formula [2] and standard Fubini theorem yield 
the following expression: 



(3.12) y(t,T) = l + Y / f y(s,T)l a3 (s,T)\ f l ai {e,T)4> H {s -9) 
~{Jo Uo 



d6 



ds. 



where y(t, T) = EexpX(t, T), X(t, T) = — £^ =1 JqI^ (*, T) 4$ for < t < T. 
By the variation of constants formula it follows that 



y(t,T)=exp(J^ e(s,T)ds 



where 



(3.13) e (i,T) = Vl ffj (f,T) [ t l (T3 (9,T) ( t )H (t-9)d8, 0<t <T. 

Therefore, we arrive at the following conclusion: condition (3.10) holds if, 
and only if, for every < T < oo 

(3.14) [ T 1 a t (y)dy = e(t,T), < t < T. 
Jo 

By differentiating expression (3.14) and making the change of variables 
x = T — t, we see that (3.10) holds if, and only if, at satisfies (3.11). □ 

In view of Corollary 3.1 we now define 

SWO == f ^ ( > • + *)<M* " 0) d9 

3=1 1 J ° 

(3.15) 



+ / 4(v)dy f 4(- + t-o)Mt-6)d0 

Jo Jo 



and we assume that the volatilities are regular enough in such a way that 

(3.16) [ T \\S H a t \\ E dt<oo 

Jo 

for every < T < oo. Indeed, it is not very restrictive to assume that the 
volatility at satisfies such integrability condition on the forward curve space 
E given in Remark 2.1. See [10], Section 5.2, for more details. 



3.2. Drift condition and quasi-martingale measure. Similar to the semi- 
martingale case, the measure P is considered as physical measure. This mo- 
tivates the following definition. 
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Definition 3.2. We say that an equivalent probability measure Q ~ P 
is a quasi-martingale measure if the discounted bond price process Z t (T) 
has Q-constant expectation, that is, for every < T < oo, 

(3.17) E Q Z t (T) = P(0,T), 0<t<T. 

Remark 3.4. This notion was already introduced in [27] in a fractional 
Black-Scholes economy under the name average risk neutral measure and 
also used in [12] in a frictionless fractional bond market in the context of 
pathwise integration. 

Remark 3.5. One should notice that the change from the actual mea- 
sure P to a quasi-martingale measure changes the mean rate of return of the 
bond but not the volatility. In fact, under a quasi-martingale measure, all 
financial assets in the bond market have the same expected rate of return, 
regardless of the "riskiness" of the bond. This is in contrast to the physical 
measure where more risky assets have a greater expected rate of return than 
less risky assets. 

We now state the main result of this section. Before this, we recall an ele- 
mentary result concerning a Girsanov theorem in the fBm setting. Without 
any loss of generality we take U = M. d . Let H be the Cameron-Martin space 
associated to the fBm. Recall that H = Image IC where 

(3.18) Kh(t) := / K{t,s)h(s)ds; h € L 2 {0,T*;R d ) 

Jo 

is an isomorphism between L 2 (0,T*;M. d ) and H. The kernel K in (3.18) is 
given by 

(3.19) K(t,s) :=c H s l l 2 ~ H j\u - s) 11 ^' 2 ^^' 2 du 

for some c# > (see, e.g., [2, 18] for more details). The next result is a 
straightforward consequence of the representation of fBm in terms of the 
standard Brownian motion. 

Lemma 3.3. Let B = . . . ,/3 rf ) be a d- dimensional fBm and let {^(t); 
< t < T*} be an M. d -valued measurable function such that J Q T \\j(t) \\ R d dt < 
oo and R(-) := J 7( s ) ds S H. Then B t := B t — J Q tr y(s) ds is a d-dimensional 
Qt* -fBm on [0, T*] such that 



(3.20) 
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where 



S^IC^R-W)^ :=exp 



{KT L R- W) T * 



2 Jo 



T* 



\\K- y R{t)\\i d dt 



and (/C 1 R • W)t* is the usual ltd stochastic integral with respect to the 
Brownian motion W associated to B. In this case, we may write 



where fi\ := (3 J t — J^ids is a Qr*-real valued independent fBm for each 
l<j<d. 



Recall that all economic activity is assumed to take place on the finite 
horizon [0,T*]. Let us fix k > which corresponds to arbitrary small propor- 
tional transaction costs in the bond market. The main result of this section 
is then the following. 

Theorem 3.1. Assume that the volatility satisfies assumptions in Lemma 
3.2 and, there exists an M. d -valued measurable function jt satisfying assump- 
tions in Lemma 3.3 in such way that 

(3.21) at7t = Suet -a t ; t> 0. 

Then there exists a quasi-martingale measure for the bond market. In addi- 
tion, the market is arbitrage free on [0,T*] with transaction costs k. 



Proof. The forward rate is the continuous mild solution of 

d 

dr t = [Ar t + a t )dt + ^a 3 t d/3 3 t , 

3=1 

under the measure P. Assuming assumptions in Lemma 3.3 and (3.21), we 
may write 

d 

dr t = (Ar t + S H o- t )dt + Yl a t d Pt 

3=1 

under the equivalent probability measure Qt* with respect to P given in 
(3.20). By (3.16) one should notice that the above equation is well defined 
under Qt* ■ By changing the measure P to Qt* in Corollary 3.1 it follows 
that for every < T < oo 

E QT ,Z t (T)=P(0,T); 0<t<T, 
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and therefore Qt* is a quasi-martingale measure. By Lemma 3.2 the two 
parameter process (t,T) \— ► X^j=i Jo -^o-j : ( s > T) dfy has QT*-fu\l support and 
therefore logZ has Qy-full support as well. By Lemma 3.1 and Proposition 
3.1 we conclude the proof. □ 

Remark 3.6. One should notice that if there exists a quasi-martingale 
measure then it should be of the form (3.20), (3.21). 

The next result gives an explicit formula for the term structure of bond 
prices in terms of a conditional expectation. By representing the fBm in 
terms of a stochastic convolution with respect to a Gaussian martingale 
M J (see, e.g., [23]) one can easily obtain an explicit formula for the term 
structure of bond prices. In the sequel, we denote 



0*(r,t):= / a j (s,t)s H - 1/2 {s-r) H - 3 / 2 ds, 



l,...,d, 



for < r < t < oo. We also write [M J ] to denote the usual quadratic variation 
of the martingale M 3 . 

Proposition 3.2. Assume that Q is a quasi-martingale measure. Then 
the bond price can be expressed by 

rT 



(3.22) P(t,T) = e 5{t ' T) E Q e 

where the kernel £(t,T) is given by 



r s {0)ds 



and G(t,T) 
u)du) 2 d[Mi] r . 



T) :=- [ T l SHa (s,T)ds-J2 fl aJ {s,T)dpi-G{t,T) 



ELiJo.fr !V(r,u)dudM} 



Proof. Formulas (2.11) and (2.10) yield the following expression: 



P(t,T) = e ?(4 ' T) E Q 



exp 



T 



r s (0)ds 



where 



at 



,T) = - fls^S^ds-Y, flaiis^dPt 

Jo Jo 



InEc 



exp 



d 

E 



T 



l aj (s,T)dPt 
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We only need to compute the above conditional expectation. By using the 
representation given in [23] it follows that 

t a j (s,t)dpi = [ d j {r,t)dMi. 



JO 

By using conditions (2.6), (2.8) and (2.9) and changing the order of integra- 
tion, we obtain 

[ T l a3 (s,T)d(3i= [ T f 9\r, U )dudMi, j = l,...,d, 

JO JO Jr 

and therefore the conditional expectation can be written as 

ft rT i _A_ r T / r T \ 2 

9 J {r,u)du) d[Mi] 



□ 



Remark 3.7. A similar pricing formula was obtained in [12], formula 
(16). One should notice a slight difference between these two formulas due to 
the fact that the drift computed in [12]; formula (15), which realizes (3.10) 
is not the correct one. See (3.13) and (3.14) for the proof of this fact. 

Remark 3.8. Besides giving an explicit expression for the term struc- 
ture of bond prices, formula (3.22) can also be used for pricing contingent 
claims in a formal way as follows. Let Q be a quasi-martingale measure and 
assume that X £ L 1 (r2,JFy, Qj<) is a claim due at time T. Then it follows 
that 



(3.23) P(t,T)E QT [X\F t }=S (t)K Q 



X 



So(T) 



a.s. Vt € [0,T], 



where the right-hand side can be formally interpreted as the price of the 
claim at time t and Qt ~ Q is defined by 

(3.24) " a -- P < T ' T > 



dQ S„(T)P(0,T) 

Formula (3.23) is very useful as long as one can compute the distribution 
of X under Qt- In the semimartingale case, this is done by means of the 
martingale property Eq[^^|^] = Zt(T) /P(0, T). In the fractional case, one 
can express this conditional expectation in terms of a prediction formula in 
the same spirit of [13]. 



We now show some examples of familiar short-rate models in the fBm 
setting as developed in this section under the quasi-martingale Q. 
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Example 1 (Ho-Lee). Let us assume that d= 1 and <7t(x) = cr, a con- 
stant, for all (t,x) € HL?.. Then in this case the model is /c-arbitrage free for 
every k > and the short-rate dynamics under Q is given by 

r t (0) = r (t) + a 2 f f '[2t - {s + 6)]<f> H {s - 0)d6 ds + afc, 
Jo Jo 

and we recognize this as the Ho-Lee model with a deterministic time- varying 
drift. 

Example 2 (Hull- White). Again assume d= 1 but now take at(x) = 
<rexp(— ax), where a and a are positive constants. Straightforward integra- 
tion imply 

r t (0)=r (t) + — /V Q (*- S ) f a [l-e- a ^]<f> H {s-e)dBda 
a Jo Jo 

+ — /'[l-e- a(i - s) ] I* e- a{t - 9) <J)H(s-9)d8ds 
a Jo Jo 

+ a f'e-^-^dps, 
Jo 

which is consistent with the Hull- White model or the Vasicek model with 
time-varying drift parameters. Furthermore, in the fBm setting this model 
is /c-arbitrage free for every k > 0. 

Remark 3.9. One can easily extend the results of this section to a 
cylindrical fBm in the HJM equation (2.2) by imposing standard growth 
conditions on the volatility in (2.6), (2.7), (2.8) and (2.9). Moreover, the 
arguments used in Lemma 3.2 are also valid in this case due to the indepen- 
dence of infinitely many fBm. 

4. Consistency for fractional HJM models. In this section, we study con- 
sistency problems related to HJM models introduced in the previous section. 
Without any loss of generality we now assume that the following dynamics 
takes place: 

d 

(4.1) dr X0 (t) = {Ar{t)+S H v(t))dt + Y,° j df3 j (t), r = x eE, 

i=i 

where Sh is given by (3.15). We study the important case of time-homogeneous 
HJM models, in the sense that the volatility a G £(U;E) does not depend 
on time t € [0, T] where < T < oo is a fixed terminal time. 

We want to stress here that the main invariance result of this section (The- 
orem 4.2) holds for any H > 1/2 when dealing with general time-dependent 
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drifts 1 1— > at £ E. We take 5//<r(t) in (4.1) only because it is the canoni- 
cal drift when calibrating HJM models from the observed forward curves. 
Moreover, it is straightforward the extension of Theorem 4.2 to the case of 
a cylindrical fBm and a Lipschitz vector field as the drift. 

Let V be the interest rate model produced by rt and let M be a param- 
eterized family of smooth forward curves (e.g., Nelson-Siegel or Svensson 
families). We recall that a pair (V,M) is consistent if all forward curves 
which may be produced by the interest rate model V are contained within 
the family A4, provided that the initial curve is in A4. There are several 
reasons, why in practice, one is interested in consistent pairs (V,M) with 
respect to HJM dynamics (see [3, 10]). In particular, the following questions 
are of great importance in calibrating interest rate models: 

Given an interest rate model V and a family of forward rate curves Ai , 
what are the necessary and sufficient conditions for consistency? Let M. 
be an exponential-polynomial family of smooth forward curves. Is there a 
nontrivial V which is consistent to Ml 

The remainder of this paper will be devoted to answer the above problems 
in the fBm setting. 

4.1. Consistent pairs (V,M). In this section, we give a fairly complete 
characterization of a given M to be consistent with respect to V. We adopt 
an abstract framework by considering M as a finite-dimensional smooth sub- 
manifold of E. The concept of invariance used in this work is the following: 



Definition 4.1. A closed set K C E is said to be invariant for the 
forward rate (r(t))o<t<T when 



Hence the pair (V,A4) is consistent if and only if M is invariant for the 
correspondent forward rate rj. The main difficulty in characterizing consis- 
tent pairs (V,M) is the obtention of the topological support for the law of 
{r(t);0 <t <T} on the space C(0,T;E) of the E-valued continuous func- 
tions on [0, T]. We recall that the topological support supp/i of a probability 
measure fj, on a Polish space X is the smallest closed set in X with total 
mass. 

Recall that the forward rate satisfies the following equation: 



F(r x °(t) eK,Vt £ [0,T]) = 1 



for every xq £ K. 



r(t) = S(t)x + S(t-s)S H o-(s)ds + Y, S(t - s)^ dft{s) 
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where {S(t),t > 0} is the right-shift semigroup acting on the Hilbert space 
E. Let us denote 

(4.2) Z(*):=£ f S{t-s)^d(3\s). 

Our main task is to characterize suppPz on the space C(0, T; E). For this 
purpose, we take advantage of the fact that fBm is a centered Gaussian 
process. The theory of Gaussian processes provides a sharp characterization 
for the support of the measure Pz ■ A direct (but lengthy) calculation shows 
that the law of Z(-) in L 2 (0,T; E) is a symmetric Gaussian measure whose 
covariance operator is given by 

^HP(t)--=[ g H (t,s)ip(s)ds, 
Jo 

where 

rsAt rsAt 

gn(t,s):= / / S(t — v)aa* S* (s — u)4>h(u — v) dudv. 
Jo Jo 

By condition (2.7) it follows that Pz is concentrated on Cq = {u £ C([0, T]; 

1 /2 

E) :u(0) = 0}. Therefore, the closure of Image in the Co-topology is the 
support of Pz- This fact would lead to a straightforward characterization of 
suppPz as long as we are able to calculate the square root of the covariance 
operator A#. In fact, a direct calculation proves to be quite difficult. More- 
over, it is very hard (maybe unlikely) to find explicitly a bounded linear 
operator A such that Ah = AA*; see Corollary B.4 in [7]. Therefore other 
nondirect techniques should be applied. 

In the sequel, we consider the Wiener space (W,H,P) of the Revalued 
fBm, where W is the space of the Revalued continuous functions / on [0, T] 
such that /(0) = 0, EI is the correspondent Cameron-Martin space defined in 
(3.18) and P is the Wiener measure on W. We have the following sufficient 
conditions for inclusions of the support of the law of an abstract Wiener 
functional V : W — ► X. See Aida, Kusuoka and Stroock [1] for the details. 

Proposition 4.1. Let V:W^ X be a measurable map, where X is a 
separable Banach space. 

(i) Let Ci : EI — > X be a measurable map, and let J n :W — ► H be a se- 
quence of random elements such that for any e > 0, 

(4.3) limP(||V-Ci° Jn\\x>e) = 0. 
Then 



suppPv CCi(EI). 
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(ii) Let C2 : H — > ^ 6e a map, and /or eac/i /iced h G H Zei T^ 1 : W — > W 
6e a sequence of measurable transformations such that ~P T h <C P for every 
n, and for any e > 0, 

(4.4) limsu P P(||V o T* - Q 2 (h)\\ x < e) > 0. 

n 

Then suppPy D (2(H)- 

The remainder of this section will be devoted to the characterization of the 
topological support of the forward rate r : f2 — > C(0, T; £7) by using conditions 
(4.3) and (4.4). Clearly, the main step is the characterization of the support 
of the probability measure P#. In the sequel, we write E to denote the 
expectation with respect to P. 

4.2. Invariance for HJM models. We now introduce a polygonal approx- 
imation for the fBm. Let us recall the Volterra representation of the fBm 

(4.5) /9(t)= /* K(t,s)dW(s); < t < T, 

Jo 

where W is the unique Wiener process that provides the integral represen- 
tation (4.5) and K(t,s) is the kernel defined in (3.19). See, for example, [26] 
for more details. 

Remark 4.1. From the above representation we notice that W is adapted 
to the filtration generated by the fBm (3 and both processes generate the 
same filtration. 

Let II = = to < t\ < ■ ■ ■ < t n = T be a partition of [0,T] where tk '■= 
and |II| := maxo<j< n - — tj) = —■ Let us consider the following 
polygonal approximations 

(4.6) (3u(t):= [ K(t,s)dW u (s)= f K(t, s)W u (s) ds, 

Jo Jo 

where 

[tj+i - tj) 

for tj < t < t»_)_i ; j = 0, 1, . . . , n — 1. 

One can check (see [18]) that for every 7 < 1 — H there exists a constant 
Ch,7 independent of II such that 

(4.7) E sup |/?n(i) - f3(t)\ < C H „\W- 

0<t<T 
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IfweW and |n| = T/n then we define wW(t) = (uj[ n) (t) , . . . , wj (t)) 
where 

J n) (t):= f K(t,s)W Ui {uj){s)ds, l<i<d. 
Jo 

Obviously u/ n ) € EI for all n > 1 and uj G W. For each /i € HI we define 

(4.8) T*u:=u + (h-uW). 

The next result is a straightforward application of a Girsanov theorem 
[18] in the fBm setting. 

Lemma 4.1. IfheW then P T h < P for all n > 1. 

It will be convenient in the sequel to make use of fractional integration 
and differentiation. For a > 0, we define the fractional integration operator 
I a and the correspondent fractional differentiation operator D a by 



1 



iaf(t) = f{aj I- (*- s r~V(*)^; 







1 d 



(4.10) D a f(t) = r - Fr ^^r / (t ~ s)- a f( S ) ds. 

1 — T(a) dt Jo 

For a comprehensive survey of the properties of these operators, see [25]. 
The most important property that we are going to use here is that I a and 
D a are each other's inverses. The following result is crucial to get (4.3) in 
Proposition 4.1. In the sequel, we write ■ (5) and • j3u) to denote the 
Paley- Wiener integrals with respect to (3 and /3n> respectively. 

Proposition 4.2. Let /3n be the polygonal approximation of the real- 
valued fBm. Ifty£L 2 (0,T;E) then 

lip E sup || (*•/?) (i)-(*-/? n )(i) lb = 0. 
l|n|Ho o<t<T 

Proof. We proceed by approximating ^ by step functions /. Assume 
that 

71-1 

f( s ) = a iX[ Sl ,s l+1 )(s) ; = so < si < ■ ■ ■ < s n = T. 

8=0 

By the semigroup property of fractional integrals and taking into account 
that D H+1 / 2 is the inverse I H+1 / 2 it follows that 

\\(f ■ m) ~ (f ■ PuMWe 
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n-1 

E 

i=0 
n-1 



a: 



((3{t i+1At )-f3(tiAt)) 



U+iM 



UAt 



Oh Puis) ds 



i+iAt) — P(U A t)) — (f3n(ti+i At) — Pu(U A *))]• 

i=0 

By the estimate (4.7) we conclude that the assertion is true for step func- 
tions. Now let us consider 4> € L 2 (0,T;E) and a sequence (f n )n>i of step 
functions which converges to ^ in L 2 (0,T; E). We have 

sup m-m)-(*-Pn)(t)\\E< sup \\(tf-(3)(t)-(f n -(3)(t)\\ E 

0<t<T 0<t<T 

+ SUp \\(fn-m)-(fn-Pn)(t)\\E 

0<t<T 

+ sup ||(/ n -/3n)(t)-(*-/8n)(t)|b 

0<t<T 

= T 1 (n)+T 2 (n,U)+T 3 (n,U). 

By the first step we only need to estimate T± and T3. Holder's inequality 
yields 

(4.11) T 3 (n,U)<\\f n -^\\ L 2 {0jT . E) \\e H Pu\\L2(o,T ] R) <00 a - s -> 

where 0jj '■ = I H ~ l l 2 oD H+1 l 2 and ||#H"/3n||L 2 (o,T;]R) is square integrable. There- 
fore we can conclude that lim^^co ET3(ra, n) = for each partition 14. It 
remains to estimate T\. For this we shall use the factorization method on 
the fractional Wiener integral. Recall the identity 



(4.12) 



7T 



sin ir a 



\Q-1, 



s - a)- a ds; a < s < t, < a < 1. 



Fix < a < 1/2 and p > l/2a. By using (4.12) and the stochastic Fubini 
theorem for fractional Wiener integrals [20], we may write 



((*-/„)■/?)(*) 



sm7ra 



7T 



\a-l 



y m (s)ds, 



where y m (a) : = / s (tf - f m )(a)(s - a)~ a d(3{a). Holder's inequality yields 

fT 



sup I 

0<t<T 



{^-fn)-P)(t)\\ 2 I<Ci / \\y m (s)\\ 2 E p ds, 



where the constant C\ depends only on p, a and T. We now choose p = l. 
The ordinary Fubini theorem and the isometry of the fractional Wiener 
integral with the reproducing kernel Hilbert space of the fBm (see [2, 9]) 
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allow us to write 

BT 2 (n)<d [ T B\\y m (s)\\ 2 E ds 
Jo 

= Cl C f f <(* - fn)(u)(s - U )- a , (tt - /„)(«)(* " 

x 4>h(u — v) dudvds. 



o Jo Jo 



By the estimate (11) in [2] we can find positive constants C2 and C3 such 
that 

E7?(n)<C a T f ||(* - /„)(«)(« - u)~ a |||dud« 

JO jO 
< C 3 || 1 I' - fn\\L 2 (0,T;E)- 

Summing up all the estimates we complete the proof of the proposition. 

□ 

We say that a closed set K in E is invariant for the evolution equation 
± y M(t) =AyM(t) +S H a(t) + aI H - l ' 2 u(t), 

(4.13) 

y(0)=x £E, 

if for each initial condition xq £ K and a control u € L 2 ([0, T]; U) we have 
y(*o,«)(t)eK for all t £ [0, T]. 

In accordance with Proposition 4.1, we are now in position to define the 
following mappings: 



(4.14) ( 1 h(t):=Y j [ S(t-s)a i 9 H h i (s)ds, heW, 

(4.15) J n (co):=J n \ w£W, 
where 6 H = I H ~ 1 / 2 o D H+1 ' 2 . 

Theorem 4.1. A closed set is invariant for the differential equation 
(4-13) if and only if it is invariant for the HJM equation (4-1)- I n particular, 



suppP^o ={y( x °> u );ueL 2 (0,T;U)}. 

Proof. We apply Proposition 4.1 to the Wiener functional Z defined 
in (4.2) with the correspondent mappings Ci, J n and T£, defined in (4.14), 
(4.15) and (4.8), respectively. Conditions (4.3) and (4.4) in Proposition 4.1 
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are direct consequences of Proposition 4.2 and Lemma 4.1. We then have 
the following characterization 

snpp Z = Ci(I H+1 / 2 {L 2 (0,T;U))) 

and therefore the law of Z is concentrated on the set of continuous functions 
of the form 

C Sit - s)aI H - 1 / 2 h(s) ds, h e L 2 ([0, T];U). 
Jo 

Now the proof follows the same lines of ([22], Proposition 1.1), and therefore 
we omit the details. □ 



4.3. Nagumo conditions and finite- dimensional invariant manifolds. In 
this section we prove the main result of this section. If M. is a C -manifold 
in E then we write T X J\4 the associated tangent space at x G M. We now 
provide Nagumo-type conditions for an HJM model to be invariant with 
respect to a given smooth manifold. 

Theorem 4.2. Let A4 be a C 1 -submanifold in E, closed as a set and 
A4 C Dom(A). Then A4 is invariant for the stochastic equation (4-1) if and 
only if 

(4.16) Ax + S H <y(t) + av G T X M 

for each x € t € [0, T] and v EU. 



Proof. Let £ be the set of [/-valued piecewise constant functions. We 
claim that a closed set K is invariant for (4.13) if and only if its mild 
solution satisfies the following condition: for each x G K and v G £ we have 
y (x,v)(t) G K for all t G [0,T]. We fix an arbitrary u G L 2 (0,T;U) and let 
us consider a sequence of step functions u n converging to u in L 2 (0,T;U). 
Then 

(4.17) \\y {x ' Un \t)-y^(t)\\ E < sup " 



0<t<T 



f S(t-s)aI H ~ 1 / 2 (u n -u)(s)ds 
Jo 



By Holder inequality we have 



(4.18) 



sup 

0<t<T 



<C 



S(t - s)aI H -^ 2 (u n - u)(s) ds 

1/2 



\<T\U r . 



u){r)\\ E dr 



where C is a positive constant which depends on T and H. Since a is 
bounded we then have inequalities (4.17) and (4.18) imply that a closed 
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set K is invariant for (4.13) if and only if y( x ' v \t) G K, t G [0,T] for all 
x £ K and all piecewise constant [/-valued function v. Thus proving our 
first claim. 

Now let M. C E be a closed C 1 -submanifold where .M C Dom(>l). By 
Theorem 4.1, Theorem 2 in [19], and the linearity of U, I H ~ l l 2 and a we 
may conclude relation (4.16) □ 

We are now in position to characterize a given finite-dimensional invariant 
submanifold. In the sequel, we write Im<7 := all. 

Corollary 4.1. Let M. be a finite- dimensional C 1 -submanifold in E 
(closed as a set). Then A4 is invariant for an HJM model given by (4-1) if 
and only if M. C Dom(A), and 

(4.19) Ax G T X M, 

(4.20) S H a{t) + lma CT X M 
for every t G [0, T] and x G M. 

Proof. Assume that A4 is invariant for an HJM model given by (4.1). 
We may apply the same arguments as in Lemma 2.3 and 2.4 in [22] to show 
that M C Dom(A) and therefore 

¥(r X0 (t) G Dom(,4),Vt G [0,T]) = 1 for every x G M. 

From Theorem 4.1 we know that every mild solution of (4.13) is also a 
strong solution such that «S#cr(0) = (j/ H ~ 1 / 2 u(0) = for every u G L 2 (0, T; U). 
By differentiating at t = and observing (4.16), we conclude relations (4.19) 
and (4.20). Conversely, let x G M, v G U and t G [0, T\. By Theorem 4.2 it 
is sufficient to check (4.16). But this is a straightforward calculation using 
the parameterizations in M. □ 

5. Nelson-Siegel family. In this section, we investigate the Nelson-Siegel 
exponential family M = {F(-,y);y G y} widely used to fit term structure of 
interest rates. The form of the curve is given by the following expression 

(5.1) F(x,y) = yi + y 2 e- y ± x + y 3 xe- y * x , x>0, 

where we restrict the parameters to the following state space y := {y = 
(yi,...,y 4 )G]R 4 |y4/0}. 

As an example, we now study simple interest rate models. By using Corol- 
lary 4.1 [in particular relations (4.19) and (4.20)], the calculations are minor 
modifications from [3] so we just give a sketch of the proofs. 
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5.1. Ho-Lee model. Let us consider the one-factor model under the Ho- 
Lee volatility structure given by a constant volatility, that is, o~t(x) = <r, for 
all (t,x) € M+. In this case, the drift restriction is given by 

(5.2) S H a(t,x) = a 2 ( 1 H (t,x) + 2 H (t)), (t,x)eR 2 +; 

where g 1 H (t,x) = 2x $Q(j) H {t-0)dO and Q 2 H {t) = t J* (j) H (t- 9) dO - J * 9<j> H {t- 
9)d9. One should notice that relation (4.19) is satisfied but because of the 
term g l H {t,x) in (5.2), relation (4.20) is not possible and therefore by Corol- 
lary 4.1 we conclude that the Ho-Lee model is not consistent with the 
Nelson-Siegel family. 

5.2. Hull-White model. Let us consider the one- factor model under the 
Hull- White volatility structure given by o~t(x) = ae~ ax where a and a are 
positive constants. In this case, the drift restriction is given by 



f (l + e - a ^)cj> H {t-9)d9 
Jo 

a(t - e UH(t-9)d9 



Sncrlt, x) = — e 

a 

(5.3) 

fe 
a Jo 

for each (t,x) € R+. By considering the full state space y, clearly the Hull- 
White model cannot be consistent with the Nelson-Siegel family. By re- 
stricting the state space to y a = {y = (y±, . . . , 2/4) I2/4 = a }) t ne curve shape 
is then given by 

F(x, y) = yi + y 2 e~ ax + y 3 xe~ ax . 

Due to the second term in (5.3), the fractional Hull- White model is not 
consistent with the Nelson-Siegel family on y a . In fact, the following result 
is not surprising in view of the previous examples. 

Proposition 5.1. There is no nontrivial fractional interest rate model 
with deterministic volatility which is consistent with the Nelson-Siegel fam- 
ily. 

Proof. By using Theorem 4.1 the proof is analogous to [3], Proposition 
7.1. □ 



Remark 5.1. One can easily find consistent families by modifying, for 
example, the Nelson-Siegel family. By using relations (4.19) and (4.20), the 
arguments are completely analogous to [3]. 
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APPENDIX: INTEGRATION FOR C R+ -VALUED PROCESSES 

In this section, we give the details of the integral introduced in Section 3 
and we keep the same notation introduced there. From Section 3 we know 
that if (J, £ V and G satisfies (HI) then J /i s dG s is well defined. In the 
sequel, tj» := iT* /2 n for i = 0, 1, . . . , 2 n ; n > 1 and < T* < oo. 

Lemma A.l. IffiEV and G satisfies assumption (HI) then 

2 n -l 



lim E 

n— >oo 



M*r(^i+iA- -G^a-) - / HsdG. 
j=o Jo 



0. 



Proof. Straightforward estimates by approximating [i by simple pro- 
cess. □ 

Let G n and G n be the upper and the lower approximations of G, respec- 
tively, along the partition (tf). Let Jq G n (s) d/i s and J^G n (s)dfj, s be the 
respective upper and lower integrals. 

Lemma A. 2. Assume that fj, € V where (H2) holds and consider G a 
stochastic process such that (HI) holds. Then: 

(a) lim^oo E sup < 4 < T , | / * G n dfx - f G n dfi\ = 0, 

(b) lim^m^oo E sup < t < T « | / * G n d\i - f G m dfx\ = 0. 

Proof. By the continuity of G and assumptions (H1)-(H2) we may 
apply dominated convergence theorem to conclude part (a). Again by conti- 
nuity, it follows that sup < s t<t* \G n (s; T) — G m (s; T) \ — > a.s. as n, m — > oo. 
Moreover, it is bounded by Csup^g^T*] 2 \G(s,T)\. Again, assumptions 
(H1)-(H2) and dominated convergence theorem allow us to conclude part 
(b). □ 

By Lemma A. 2 we shall define 

/ G s dfi s := lim / G n (s)dfi s = lim / G n (s)dfi s . 
The next result is a straightforward integration by part formula. 
Proposition A.l. Assume that assumptions (HI) and (H2) hold. Then 

(A.l) / G s d^L s + [ fi s dG s = G T *HT* -G fj, . 

Jo Jo 
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Proof. By writing a telescoping sum we have 



2 n -l 



i=0 



= Gt*^t* - Go fj,o 






a.s. for all n > 1. By Lemmas A.l and A. 2 we only need to show that the 
left side goes to zero as n — > oo. But this is an immediate consequence of 
hypotheses (H1)-(H2) together with the continuity of G. □ 
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